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A GENERAL STAGNATION-POINT CONVECTIVE -HEATING EQUATION 
FOR ARBITRARY GAS MIXTURES 

By Kenneth Sutton and Randolph A. Graves, Jr. 

Langley Research Center 

SUMMARY 

The stagnation -point convective heat transfer to an axisymmetric blunt body for 
arbitrary gases in chemical equilibrium was investigated. The gases considered were 
base gases of nitrogen, oxygen, hydrogen, helium, neon, argon, carbon dioxide, ammonia, 
and methane and 22 gas mixtures composed of the base gases. In the study, enthalpies 
ranged from 2.3 to 116.2 MJ/kg, pressures ranged from 0.001 to 100 atmospheres, and 
the wall temperatures were 300 and 1111 K. 

A general equation for the stagnation- point convective heat transfer in base gases 
and gas mixtures was derived and is a function of the mass fraction, the molecular weight, 
and a transport parameter of the base gases. The relation compares well with present 
boundary -layer computer results and with other analytical and experimental results. In 
addition, the analysis verified that the convective heat transfer in gas mixtures can be 
determined from a summation relation involving the heat-transfer coefficients of the base 
gases. The basic technique developed for the prediction of stagnation-point convective 
heating to an axisymmetric blunt body could be applied to other heat-transfer problems. 

INTRODUCTION 

Convective heating during planetary entry depends on the atmospheric gas composi- 
tion surrounding the planet. For the planets presently considered for entry, the atmo- 
spheres are primarily binary mixtures of nitrogen and oxygen (air) for Earth; carbon 
dioxide and nitrogen for Mars and Venus; carbon dioxide and argon for Mercury; and 
helium and hydrogen for Jupiter. (See refs. 1 to 7.) Trace amounts of neon, methane, 
and ammonia have also been considered for several of the planets. Except for Earth, the 
actual composition of the atmospheres are not accurately known. Consequently, numerous 
calculations will be required to determine the convective heating during planetary entry 
unless a general relation can be determined. Generally, existing analytical (refs. 8 to 16) 
and experimental (refs. 17 to 27) studies have emphasized only particular gases or gas 
mixtures and a solution for arbitrary gases was not obtained. In references 14, 16, 


and 28, general relations for convective heating are presented but the analyses were 
developed for a limited number of gases and gas mixtures. 

The purpose of the present study is to develop a general relation for calculating the 
convective heating to the stagnation point of a blunt axisymmetric body for gases or gas 
mixtures in chemical equilibrium, especially those gases encountered during high-velocity 
planetary entry. In the analytical study, the approach is to derive a simple approximate 
relation from the basic theoretical definition of heat transfer in gases. Also, the convec- 
tive heating for numerous gases is calculated by an existing computer code (computer 
study) to obtain results for use in developing the approximate general relation and as a 
basis of comparison with the approximate relation. For the computer study, stagnation 
enthalpies ranged from 2,3 to 116.2 MJ/kg, stagnation pressures ranged from 0.001 to 
100 atmospheres (1 atmosphere equals 101.325 kN/m 2 ), and wall temperatures were 300 
and 1111 K. The gases considered were base gases of nitrogen, oxygen, hydrogen, helium, 
neon, argon, carbon dioxide, ammonia, and methane and 22 mixtures composed of the base 
gases. Radiative heating was neglected in the present study. 

SYMBOLS 


c mass fraction 

Cp specific heat at constant pressure, j/kg-K 

D reference diffusion coefficient, see equation (6), m 2 /s 

D e ff effective diffusion coefficient defined by equation (16) 

Dij binary diffusion coefficient, m 2 /s 

E energy factor, see equation (C2) 

F diffusion factor, see equation (6) 

h enthalpy, MJ/kg 

h° standard enthalpy of formation, MJ/kg - 

j species mass flux, kg/m2-s 

K heat-transfer coefficient defined by equation (33), kg/s -m 3 / 2 - atm 1/2 
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k Boltzmann constant, 1.380622 x 10"^3 J/K 

M molecular weight 

NLe Lewis number defined by equation (15) 

Np r Prandtl number defined by equation (13) 

Ng c Schmidt number defined by equation (14) 

p pressure, atm (1 atmosphere equals 101.325 kN/m2) 

q convective heating rate, MW/m2 

R equivalent body radius, m 

tfL universal gas constant, 8.31434 X 10^ J/kmole-K 

r body coordinate, m 

s distance along body surface from stagnation point, m 

T temperature, K 

T* nondimensional temperature defined by equation (B2) 

u velocity component parallel to body surface, m/s 

y normal distance through boundary layer from wall surface, m 

Mcj 

a constant, see equations (36) and (38), 2.2621 x 10 

/3 velocity gradient defined by equation (23), s~* 





y transport parameter defined by equation (37), A - ^ K”^‘^ 

e maximum energy of attraction, J 

V transformed coordinate defined by equation (20) 

X thermal conductivity, J/m-s-K 

M viscosity, N-s/m^ 

| transformed coordinate defined by equation (21) 

p density, kg/W* 

a molecular collision diameter, A (1 angstrom equals 10 “ 10 meter) 

(2 2 )* 

’ ' reduced collision integral for viscosity 

Subscripts: 

e external edge of boundary layer 

i ith species or component 

j jth species 

max maximum 

o initial composition 

s stagnation condition at boundary -layer edge 

w wall value 

00 free-stream value 
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DESCRIPTION OF COMPUTER PROGRAM 


As part of the present study, the convective heating of the base, gases and gas mix- 
tures was calculated by an existing computer code which is a numerical solution of the 
multicomponent boundary-layer equations for gases in chemical equilibrium. Results 
from these computations were used in developing the approximate general relation and as 
a basis of comparison with the approximate relation. Descriptions of the computer code 
and method of solution of the boundary -layer equations are presented in references 29 
to 32. The routines for the thermodynamic and transport properties are so constructed 
that any arbitrary gas mixture can be considered. Hereinafter, these solutions are 
referred to as computer results. 

The thermodynamic properties are calculated by the method described in refer- 
ence 33, and, for the present study, the necessary inputs were obtained from references 34 
to 37. The inputs are for a common basis, and a comparison of this method with other 
methods is given in reference 33. 

The mass diffusion is calculated by using the bifurcation approximation as presented 
in reference 38. The accuracy of the method depends on the bifurcation diffusion factors, 
and, for the present study, the factors are calculated by the method given in reference 39. 
For a few species, the basic data required were not available, so the diffusion factors 
were taken from the correlation of reference 38. Comparisons of the bifurcation diffusion 
approximation with exact multicomponent diffusion are given in references 30, 38, and 39. 

The viscosity and thermal conductivity values are calculated by approximate rela- 
tions of the Sutherland-Wassiljewa type as presented in reference 38. The method 
requires the previously mentioned bifurcation diffusion factors and also the self-diffusion 
factors calculated by the exact method presented in reference 34. Comparisons of the 
transport properties calculated by this method and other methods are given in 
reference 34. 

The computer program provided solutions for arbitrary gas mixtures and thus gave 
a consistent base for developing and evaluating the approximate relation. In the present 
study, the computer program was used to obtain numerical solutions for the stagnation- 
point convective heating for nine base gases and 22 gas mixtures. They are as follows - 

Base gases: 


He CO 2 

Ne NH 3 

ch 4 


n 2 

02 

H 2 


Ar 
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Gas mixtures (composition by mass fraction): 


0.2320 02-0.7680 N 2 (air) 

0.6714CO2 -0.3286 Ar 

0.1345 CO2-0.8655N2 

0.5000 002 -0. 5000 N2 

0.8500 CO2-0. 1500 N2 

0.7500 02-0.2500 N 2 

0.5000 02-0.5000 N2 

0.1500H2-0.8500 He 

0.3500 H 2 -0.6500 He 

0.6500 H2-0.3500 He 

0.5000 H2 -0. 5000 Ar 

0.1500 N2-0.8500H2 

0.5000 N2-0.5000 H2 

0.2000 CO2-O.8OOO H2 

0.4000 CO2-O.6OOO H2 

0.6000 CO2 -0.4000 H2 

0.8000 CO2-O.2OOOH2 

0.3626 CO2-0.3077N2-0.3297 Ar 

0.1339 C0 2 -0.8525N2-0.0136Ar 

0.3000 Ne-0.3000 Ar-0.4000 He 

0.3500 H 2 -0.3500 C02-0.3000 N 2 

0.352 H 2 -0.423 He-0. 176 Ne-0.024 CH4-O.025 NH3 


The 81 chemical species considered in the computer study are as follows: 


Ar 

C 3 

NO 

co 2 ' 

C 

c 3°2 

no 2 

He + 

CH 

c 4 

no 3 

He ++ 

CHNO 

c 4 n 2 

n 2 

H + 

ch 2 

C 5 

n 2 h 2 

H" 

ch 2 o 

He 

n 2 o 

HCO + 

ch 3 

H 

o 

H 3 o + 

ch 4 

HCN 

OH 

Ne + 

c 2 h 

HCO 

02 

Ne ++ 

c 2 h 2 

HNO 

°3 

N + 

c 2 h 4 

hno 2 

e“ 

N ++ 

c 2 h 4 o 

hno 3 

Ar + 

NO + 

CN 

ho 2 

Ar ++ 

no 2 " 

CNN 

h 2 

C + 

n 2 + 

cn 2 

h 2 o 

c ++ 

o + 

CO 

Ne 

C~ 

o' 

co 2 

N 

CN + 

OH + 

C 2 

NH 

CN" 

OH" 

c 2 n 

nh 2 

CO + 

o 2 + 

c 2 n 2 

nh 3 

c 2 ' 

o 2 _ 

c 2 o 





For each base gas or gas mixture, the heating was calculated for a combination of condi- 
tions of — 

Stagnation enthalpy: 2.3, 4.6, 9.3, 13.9, 18.6, 23.2, 27.9, 32.5, 37.2, 46.5, 

58.1, 69.7, 81.4, 93.0, 104.6, and 116.2 MJ/kg 

Stagnation pressure: 0.001, 0.01, 0.1, 1.0, 10.0, and 100.0 atm 

Wall temperature: 300 and 1111 K 
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ANALYSIS 


A simple general expression for the convective heating to the stagnation point of a 
blunt axisymmetric body for high-enthalpy gases in chemical equilibrium was derived 
from the basic theoretical definition of heat transfer in gases. To meet the criterion of 
simplicity, it was necessary to approximate several of the more complex expressions that 
occurred in the theoretical development with more tractable expressions. 

From reference 40, the convective heating to a wall from a gas is 



The first term on the right-hand side of the equation is the heating due to conduction and 
the second term is the heating due to diffusion. The mass diffusion due to thermal gradi- 
ents is small in comparison with mass diffusion due to concentration gradients (ref. 40) 
and has been neglected. The specific heat and enthalpy of a mixture are defined as 

c p = X C i c P,i 

h = ^ c i h i ( 3 ) 

where 

hi = j* c Pji dT + h? (4) 

With these definitions, equation (1) becomes 



The species mass flux is based on a bifurcation approximation to binary diffusion 
coefficients. This approximation has been used in references 29, 30, 38, and 39 and was 
shown to be in good agreement with exact multicomponent diffusion. The present form of 
the bifurcation approximation is taken from reference 38. The approximation is 
expressed as 


D-- * D 

Dl J Tp- 


( 6 ) 
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where D is a reference diffusion coefficient and Fj and Fj are diffusion factors for 
the ith and jth species, respectively. 

As given in reference 38, the species mass flux for the bifurcation approximation 
can be expressed as 


pD/3 2 8z i 

For unequal diffusion, the quantity zj lies between a mass and a mole fraction and is 
defined by 


z i = 


MCj 


( 8 ) 


and /3j and ^ are system quantities defined by 


^2 = M 2 


c i F i 


From equation (8), 


M_^i = 0 8z i ! „ c i 8M 

Fj By 2 9 y 1 9 y Fj 9 y 

V -V" ^ 

Neglected 


(9) 

( 10 ) 


( 11 ) 


The last two terms on the right-hand side of equation (11) are neglected in the present 
analysis. The neglected terms involve variations in system quantities and the variations 
of these quantities should be much less than the variations of individual species. Thus, 
in the present analysis, it is assumed that 

9&i F^o 9z, 

— i (12) 

dy M 9y 


With the following nondimensional parameters — 

Cpp. 

Prandtl number: Np r = —— 

A 


(13) 
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Schmidt number: 


(14) 


n ct - a 
SC pDeff 


Lewis number: 


pCpD e ff Npr 


where 


D eff = 


and equations (5), (7), and (12), the convective heating can be expressed as 

q = - — y Fihi — + N Le Y ^ ^ 

W N Pr,w\fy M L 11 8y Le ^ dy) 


In the present study, an analogy between energy and mass transfer is used in the 


form 


i _ 1 i,e ~ ^i/w ah 

y ^ 1/3 h g - h w 9y 


This analogy is based on the analyses presented in references 30, 38, 40, 41, and 42. By 
using equation (18), equation (17) becomes 


. _ /ah\ 1 02, w ^ F i h i,w( z i,e - z i,w) / \2/3 ^ h i,w( z i,e - z i,w) 

qw " N Pr w \9yj " / \l/3 h s - h w ' Le,W ^ h s - h w 

w M w (N Le>w j 


The present equation is transformed into the usual boundary-layer coordinate sys- 
tem for a blunt axisymmetric body by the Lees-Dorodnitsyn transformation (ref. 43): 


pc* 
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( 21 ) 


* = j* 0 Pe^e u e r2ds 


For a stagnation-point solution 

P e M e u e 2r2 


I =■ 


4/3 


( 22 ) 


where the velocity gradient @ is 


P - 


du e 

ds 


(23) 


Thus, the transformation from the y coordinate to the tj coordinate for the stagnation- 
point solution of a blunt axisymmetric body is 


V 


|P e P e J 0 


(24) 


and applying this transformation, equation (19) becomes for the stagnation point 


• _ P 1//2 foh 

q w ~ 


N Pr(PePe) 1/2 ^ /w 


^2,w ^ *ihi,w( z i,e " z i,w) 

[ w( N Le,w) 


,1/3 


hs _ 


, / KT \ 2 / 3 /^ hi M z i,e - z i,w) 
\ Le ’ W ' h s - h w 


(25) 


At the stagnation point, the velocity gradient for Newtonian flow is 




2p, 


e/j Pqq ' 


(26) 


For 


p„ » p 

^e F oo 
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the velocity gradient can be approximated by 



With the use of the ideal gas equation 


( 27 ) 


P = 


T 

M 


and the boundary -layer approximation at the stagnation point 


(28) 


Pe = Pw = p s 


(29) 


equation (25) becomes 



8 1//4 /M e 
^Pr,w\^ T e 


1/4 


P W P 


W^W 


Pe^e 


1 / 2(817 



1 - 




(30) 


An expression is now needed which includes the enthalpy gradient at the wall, and it 
was necessary to use a correlation of the results from the computer study. A discussion 
of the enthalpy gradient at the wall is given in appendix A. A correlation which is reason- 
ably valid for all the base gases and gas mixtures is presented (eq. (A3)) as follows: 



Special note should be taken of the parameter p 0 e , which is the viscosity evaluated at 
the temperature at the outer edge of the boundary layer (T e ) for a mixture based on the 
initial composition of the cold mixture or the ambient composition of the atmosphere sur- 
rounding a planet. For example, in determining the heat transfer in a gas mixture of 
0.60 N 2 -O .20 Ar-0.20 CO 2 , the viscosity e is based on this composition rather than 
the chemical equilibrium composition at the outer edge of the boundary layer. 

From the enthalpy-gradient correlation of equation (31), the heat transfer from 
equation (30) becomes 
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The stagnation-point convective heating to a blunt axisymmetric body is often expressed 
in terms of a heat-transfer coefficient. This coefficient is considered to be essentially 
constant over a range of conditions for a particular gas or gas mixture. As presented in 
references 16, 25, 28, and 44, the heat-transfer coefficient K is defined as 



When equations (32) and (33) are compared, it is seen that the heat-transfer coefficient 
can be written 


K = 


i/a/ \ 1/4/ \V8 

0.58(8) \ Mo\ , ,1/2 

, xO.eUTe) \M e y ^°> e ) 

(^Pr,wj X ! \ ' 


1 - 




(34) 


The viscosity is a dominant factor for the convective heat transfer of gases; thus, 
in order to find a general expression for all gases, the viscosity must be expressed in 
terms of basic parameters. For gas mixtures, the exact calculation of the viscosity is 
too complex to be used in any simple expression for the heat transfer. For example, see 
references 45 to 48. In the present analysis, a simple summation is proposed for mix- 
ture viscosity which agrees adequately with experimental results and theoretical models. 
The viscosity of the pure components is based on the analysis presented in reference 47. 
A discussion of the present viscosity relations is presented in appendix B. 

The relation proposed for the mixture viscosity is (eq. (B6)) as follows: 


i=M*/ 2 y 

M Li 


Mi 


1/2 


H 


( 35 ) 


13 



For a pure gas component 


Mi = < *V i M i 1/2 T 0, 


65 


where the transport parameter y ^ is 





and the constant a is 

a = 2.2621 x 1CT 6 


(36) 


(37) 


(38) 


With the use of the above viscosity relations, equation (34) becomes 



As previously mentioned, the heat-transfer coefficient is used as a constant over a 
range of conditions. In the present analysis, the Lewis number augmentation factor 
(appendix C), 



is assumed to be approximately equal to 1. This approximation is based on a comparison 
with the equivalent expression for equal binary diffusion coefficients. (See refs. 8 and 49.) 
Furthermore, the following approximation is made 


'M e \ 1//8 0.075 _ 

TO e 


( 40 ) 
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The terms in equation (40) have small exponents and will have a small variation over a 
range of conditions. Also, M e will decrease with increasing T e due to dissociation 
and ionization. The value used in the approximation was determined from the computer 
results. The effect of the parameters in equation (40) on the heat-transfer coefficient 
will be discussed further in the next section. 


With the above approximations and evaluation of the constants, equation (39) for the 
heat-transfer coefficient can be expressed as 



The Prandtl number at the wall varied from 0.39 to 0.71 for the various base gases and 
gas mixtures considered in the present study. However, the Prandtl number at the wall 
only varied from 0.67 to 0.71 for those base gases and gas mixtures which are presently 
considered for planetary entry. For Np r w = 0.69, 


K = 0.1106 




M o,iY 0>i J 


(42) 


The present analysis has shown that the approximate stagnation-point convective heating 
for gases has a general relation involving the mass fraction (c 0> i), molecular weight 

(M 0) i), and transport parameter (y Q ^ of the base gas or the base gases of a mixture. 

The necessary values for calculating y . from equation (37) are given in reference 50 

for numerous gases. Listed in table I are the values for the base gases considered in the 
present study. The approximate convective heating to the stagnation point of a blunt axi- 
symmetric body can be calculated by the use of equation (33) and either equation (41) 
or (42) for base gases and gas mixtures in chemical equilibrium. 

The heat-transfer coefficient for gas mixtures can also be calculated if the heat- 
transfer coefficients (KjJ of the base gases are known. By manipulating equation (41), it 
can be shown that 
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I 


and for N pr w - Np rwi 


_1_= V 

K 2 '^ 2 


( 44 ) 


The present analysis is for the stagnation-point convective heating to an axi sym- 
metric blunt body and the developed relation is only applicable to this problem. However, 
the basic technique used in the present development could be applied to other heat-transfer 
problems as well. 

RESULTS AND DISCUSSION 


The present analysis has shown that the convective heating for gases can be 
expressed by a general relation involving the mass fraction molecular weight 

(M 0 , and transport parameter (y Q These parameters are known for most gases. 

The results of this general expression are compared with the results of the numerical 
solution to the boundary-layer equations from the present computer study and with the 
results of analytical and experimental studies presented in the literature. 

Results from the present computer study are presented in figure 1 in the form of 

as a function of h s - h w . The linear curve shown in figure 1 for a particular gas 

or gas mixture is based on a single heat-transfer coefficient for the gas as determined by 
a least-square fit of the values of K from each data point. The heat-transfer coefficient 
for each gas could be correlated with good accuracy by a single coefficient for a range of 
enthalpies, pressures, and wall temperatures. The heat-transfer coefficients for the 
gases and the average error and maximum error for the correlation are listed in table n. 

The present analysis showed that the heat-transfer coefficient for gases can be 
expressed as (see eq. (41)) 



K = 


0.0885 


( N Pr,w) 



- 1/2 


“04 


M °,i>o,i 


and for N pr w = 0.69 (see eq. (42)) 
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The analytical expressions (eqs. (41) and (42)) for the heat-transfer coefficient are com- 
pared with the present computer results in figure 2. As shown by the comparison, the 
general expression is in good agreement with the present computer results. Furthermore, 
the data presented in figure 2(b) show that equation (41), which includes the wall Prandtl 
number, is in better agreement with the computer results. Heat-transfer coefficients 
were determined for various gases from experimental and analytical studies presented 
in the literature; and, as shown by the comparison in figure 3, there is good agreement 
between the general relation and the results from other studies (refs. 13 to 22, 24 to 26, 
and 28). In references 14 and 15 the convective heating for several base gases and gas 
mixtures was calculated by the same method. The heat-transfer coefficients determined 
from these references are also compared separately with the general relation in figure 4. 
As shown by the data, the heat-transfer coefficients from references 14 and 15 are linear 

/ x-l/2 


functions of the parameter 



“°,i 


M n . 

OA'O,! j 


; however, the slope is slightly greater than 


given by the present results. 


The previous discussions have been concerned with the validity of the final, overall 
expression from the present analysis. The trend of the variations of the heat-transfer 
coefficient can also be explained from the analysis. From the computer results presented 


in figure 1, the reduced-heating parameter 



for some gases (for example, ammo- 


nia) appears to have an "apparent” wall temperature effect. This apparent temperature 
effect is actually a wall Prandtl number effect. The analytical analysis, see equation (41), 
shows such an effect and it can be substantiated from analysis of the computer results. 

For the gases which showed a wall temperature effect, there was a variation in the wall 
Prandtl number for the different wall temperature values. 


The computer results in figure 1 also show a decrease in reduced-heating parameter 
with decreased pressure for a given enthalpy. This effect can be explained by examining 
equation (40). The value of 1.81 was based on an average over a range of conditions. For 
a given enthalpy, the degree of dissociation and ionization will be greatest at the lower 
pressures, and the molecular weight M e will be decreased. Also, this effect will cause 



the temperature T e to be lower at the lower pressures. Without the approximation of 
equation (40), equation (42) can be expressed as follows for Np r >w = 0.69: 


/ \!/8 / \ 

K = 0.0632p£) T e °' 075 | ) — 

W ^ M o,i y o4 


■ 1/2 


\ 


(45) 


The trend of M e and T e with pressure, from the computer results, is shown in fig- 
ure 5. Also, in figure 5, the heat-transfer coefficient as calculated by equation (45), with 
the value of M e and T e taken from the computer solution, is compared with the heat- 
transfer coefficient from the computer results. As shown by the comparison, the present 
analysis predicts correctly the effect of pressure. 


Equation (45) also shows the trend of the heat-transfer coefficient with enthalpy. 
The heat-transfer coefficients from the computer study and from equation (45) are com- 


pared in figure 6 for several cases. The term 


(Me \ 1/8 
M, 1 


oj 


T 0.075 gjygg ^ predominant 


trend of the coefficient with enthalpy and, as shown in figure 6, good agreement generally 
exists between the present analysis and the computer results. The present analysis gives 
good final results and also gives the correct trend of intermediate steps. 


The present results of the analytical and computer studies were compared with 
three other methods for calculating the convective heating of gases. These methods are 
presented in references 14, 16, and 28 and are different from each other and from the 
present analysis. These methods were based on a smaller number of gases than the 
present study. 


A simple method for determining the convective, stagnation-point heating of gas 
mixtures is presented in reference 28. In this method the heat-transfer coefficient of 
gas mixtures is calculated by knowing the coefficient of the base gases by the relation 


1 _ V C o,i 
K L Kj 


(46) 


The relation was derived by empirical means and was based on mixtures of CO 2 , N 2 , and 
Ar. Similar expressions (see eqs. (43) and (44)) were derived for a wider range of gas 
mixtures in the present analysis: 
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A comparison between equations (43), (44), and (46) and the present computer results is 
given in figure 7 for some binary gas mixtures and in table m for all the gas mixtures 
considered in the present study. In comparing the summation relations, the heat-transfer 
coefficients for the base gases are taken from the computer results. The summation 
relations from the present studies, equations (43) and (44), have good agreement with the 
computer results for their respective conditions. Equation (43) has the best agreement 
for all the gases; however, equations (44) is equally valid for the gases when the Prandtl 
numbers at the wall are approximately equal for both the mixtures and the component 
base gases. The empirical method of reference 28, equation (46), has good agreement 
with the computer results for the mixtures where the Prandtl numbers are approximately 
equal for the mixtures and component base gases. However, this method is not as accu- 
rate as equation (43) when the Prandtl numbers for the mixtures are different from those 
for the component base gases. The method of reference 28 was based on empirical 
means, whereas the present relations were derived from analytical considerations. From 
equations (44) and (46) and analysis of the computer results, it can be seen that differenti- 
ation between a first and second power for the heat-transfer coefficient would be difficult 
by empirical means. 

The method presented in reference 14 and later extended in reference 15 is based 
on knowing the low-temperature properties of the gas. In using the method, a correlation 
would have to be made for each particular gas, and, for air, N2, CC>2, Ar, H 2 , and several 
N 2 -CC> 2 -Ar mixtures, the necessary correlative parameter are listed in reference 15. 
Basically the method of reference 14 and the present theoretical analysis are in agree- 
ment concerning low-temperature properties. In the present analysis, the heat-transfer 
coefficient is a function of the parameters Np r w , y Q> M 0 , and c Q . These parameters 

are low-temperature properties of the gases. The low-temperature properties give the 
major component of the convective heating, and the boundary-layer-edge properties (high- 
temperature properties) give minor variations to the heating. As previously shown in 
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figure 4, heat-transfer coefficients can be determined from the results of references 14 

/ c o1 V I/2 

and 15 and were shown to be linear functions of the parameter I V — . This 

l Zy Mo,i^o,i/ 

is in agreement with the present analysis. The present summation relation of equa- 
tion (44) was applied to the data of references 14 and 15. As shown in table IV, the heat- 
transfer coefficients determined by the present relation is in good agreement with those 
determined with the use of references 14 and 15. The Prandtl number at the wall for 
gases considered for planetary entry is approximately 0.69; therefore, the present method 
is easier to use than the method presented in reference 14. 

The results of reference 16 are presented as a correlation of the heat-transfer 
coefficient with molecular weight in the form 

K = 0.0323 + 0.00233M o (47) 


or 


K = 0.0323 + 0.00233 




(48) 


A comparison of the results presented in reference 16 and the present computer results 
is shown in figure 8. The molecular weight correlation of reference 16 was based on a 
limited number of gases and for these gases the correlation is valid. However, as shown 
by the present results, this method is not valid for arbitrary gases. 

CONCLUSIONS 


The present analysis has resulted in the development of an approximate method for 
calculating the convective stagnation-point heating to an axisymmetric blunt body for gases 
in chemical equilibrium. The heating is expressed by a general relation of the mass frac- 
tions, molecular weights, and transport parameters of the base gases. Comparisons with 
the more exact computer study and with other analytical and experimental studies have 
shown the method to be valid. Also, it has been shown that a summation relation derived 
from the analysis is valid for calculating the heating of gas mixtures from the data of the 
base gases comprising the mixtures. For the gases presently considered for planetary 
entry, the approximate convective heating can be easily calculated. Since the analysis 
was based on a wide range of base gases and gas mixtures, the results should be valid 
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for most gases. The basic technique developed for the prediction of stagnation-point 
heating to an axisymmetric blunt body could be applied to other heat -transfer problems 
as well. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 8, 1971. 
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APPENDIX A 


ENTHALPY GRADIENT AT THE WALL 


The results presented in references 51 and 52 have shown that the enthalpy gradient 
at the wall for the stagnation point of a blunt axisymmetric body in high-enthalpy air can 
be expressed by 


1 fdh \ _ / Pe^e ' 
h s " h w \ v w v W^W; 


where a and b are given as follows: 


(Al) 


b Determined 

a from reference - 

0.51 0.975 51 

.49 .938 52 


Also, from analysis of the heat -transfer relations presented in references 8 and 53, 



where a, b, and c are given as follows: 

Determined 

c from reference - 

0.4 8 

.4 53 

The interest in the present study was whether a single expression similar to equation (A2) 
could be determined and whether it would be reasonably valid for all gases. If the expo- 
nent for the density -viscosity product was 1, then the heat-transfer equation (eq. (30)) 
could be made simpler by the elimination of terms. A similar approach was attempted 
in reference 14, but a single correlation could not be determined for air, nitrogen, hydro- 
gen, carbon dioxide, and argon. Preliminary work in the present study showed that such 
a correlation was feasible and additional studies were performed. 


a b 

0.54 0.900 

.54 .980 
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APPENDIX A — Continued 


The data used to correlate the enthalpy gradient were taken from the results of the 
present computer study. The initial correlation is shown in figure A1 in the format of 
equation (A2) with b = 1 for the primary base gases and air. The data of the various 
gases are in good agreement with each other; however, a straight line could not be fitted 
through the zero point and there is a slight curvature to the data. Also, the relation was 



Air 

Nitrogen 

Oxygen 

Hydrogen 

Helium 

Argon 

Carbon dioxide 


Figure Al.- Initial correlation of the enthalpy gradient. 


not in good agreement for all gases considered in the present study. Variations to the 
basic form of equation (A2) were tried and the best correlation obtained to the data was 



As shown in figure A2(a), the final correlation, equation (A3), has good agreement for the 
gas mixtures which are presently considered in planetary entry and for the base gases 
which comprise these mixtures. Representative data for all the gases considered in the 
present study are also compared with equation (A3) in figure A2(b). The greatest devia- 
tion from equation (A3) is for the gases with a significant variation of the molecular weight 
at the wall for a range of conditions, especially with stagnation enthalpy. The variation of 
molecular weight is affected by chemical reaction and diffusion. 
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APPENDIX A - Continued 



o Air 
□ Nitrogen 
O Oxygen 
a Hydrogen 
k Helium 
& Argon 
q Ca'rbon dioxide 
o 0.35 H 2 - 0.65 He 
0 0.65 H 2 - 0.35 He 
A 0.50 C02- 0.50 N 2 
© 0.85 C0 2 - 0.15 N 2 
e 0.6714 C0 2 - 0.3286 Ar 
<$> 0.3626 C0 2 - 0.3077 N 2 - 0.3297Ar 
A 0.1339 C0 2 - 0.8525 N 2 - 0.0136 Ar 


F 



o Air 
n Nitrogen 
O Oxygen 
a Hydrogen 
t\ Helium 
b Argon 
□ Carbon dioxide 
o 0.35 H2- 0.65 He 
0 0.65 H 2 - 0.35 He 
a 0.50 C0 2 - 0.50 N2 
© 0.85 C0 2 - 0.15 N 2 
0 0.6714 C0 2 - 0.3286 Ar 
<3> 0.3626 C 02 - 0.3077 N 2 - 0.3297 Ar 
A 0.1339 C 02 - 0.8525 N 2 ' 0.0136 Ar 
k Neon 
k Methane 
® Ammonia 
® 0.50 H 2 - 0.50 Ar 
$ 0.30 Ne - 0.30 Ar - 0.40 He 
a 0.60 C0 2 - 0.40 H 2 


(b) All gases considered in present study. 

Figure A2.- Final correlation of the enthalpy gradient. 


J 
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APPENDIX A - Concluded 


The correlation of equation (A3) for the enthalpy gradient at the wall was based on 
boundary -layer results for a wide range of flow conditions and a wide range of base gases 
and gas mixtures. The base gases included monatomic, diatomic, and polyatomic gases. 
The gas mixtures were composed of base gases of only high molecular weights, of only 
low molecular weights, of high and low molecular weights, of basically inert gases, and 
of reactive gases. Since the correlation was in agreement for a wide range of base gases 
and gas mixtures, the correlation should be equally valid for most gases. 



APPENDIX B 


VISCOSITY RELATIONS 


The viscosity equation for a pure component from reference 47, is 


Mi = 2.6693 X 10" 6 



(Bl) 


where 

Mj viscosity, N-s/m^ 

T temperature, K 

Mi molecular weight 

<7i collision diameter, A 

(2 2 )* 

O' 5 ' reduced collision integral for viscosity 

The reduced collision integral is based on the Lennard-Jones (6-12) potential 

and is a slowly varying function of the nondimensional temperature T* defined as 

T* » £ (B2) 

where e is the maximum energy of attraction and k is the Boltzmann constant. The 
values of as a function of T* are presented in reference 47 and are shown in 

figure Bl. For the species and temperatures of interest in the present study, the reduced 
collision integral can be expressed as 

fi (2,2)* = 1.18(T*r 0 * 15 (B3) 

By combining equations (Bl), (B2), and (B3), the viscosity for a pure species can be 
expressed as 


M t = (2.2621 X 10-6)y.M i 1//2 T 0 - 65 


(B4) 
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APPENDIX B - Continued 



Figure Bl.- The reduced collision integral for viscosity as a function 
of the nondimensional temperature. 

where the transport parameter is given as 


d = 


•ft) 


0.15 


(B5) 


A tabulation of the collision diameter a- and the maximum energy of attraction e^k 
for numerous species are given in reference 50. 

A simple summation expression is proposed for the mixture viscosity. This 
expression is given by 


i = ul/2 l 


Mi l/2 Mi 


(B6) 


The exact treatment of mixture viscosity involves more complicated expressions. (For 
example, see refs. 45 to 48.) The simple relation given by equation (B6) has been com- 
pared with experimental results and with more exact calculations. Some results of the 
comparison are shown in table BI and figure B2. (Also see refs. 54 to 60.) As shown by 
these results, equation (B6) is in good agreement with experimental and calculated 
values. 


27 



APPENDIX B - Continued 


TABLE BI.- COMPARISON OF MIXTURE VISCOSITIES 


Gas composition 
(composition by mole fraction) 

Temp., 

K 



p, N-s/m^ 

Reference 

From 

reference 

Computed (eq. (B6)) 
with from reference 

0.1754 He - 0.5556 Ne - 0.2670 Ar 

293 

54 

27.34 X 10 -6 

27.05 x 10 ' 6 

0.5175 He - 0.3210 Ne - 0.1615 Ar 

293 

54 

25.86 

26.42 

0.3594 He - 0.3193 Ne - 0.3213 Ar 

473 

46 

36.00 

36.65 

0.3333 H 2 - 0.3333 Ne - 0.3333 C0 2 

298 

46 

19.02 

18.72 

0.2500 H 2 - 0.2500 Ne - 0.2500 C0 2 - 0.2500 CC1 2 F 2 

298 

46 

15.93 

16.69 

0.2500 H 2 - 0.2500 N 2 - 0.2500 C0 2 - 0.2500 CC1 2 F 2 

298 

46 

14.67 

15.55 

0.1754 He - 0.5576 Ne - 0.2670 Ar 

473 

55 and 57 

37.90 

38.08 

0.3594 He - 0.3193 Ne - 0.3213 Ar 

473 

55 and 57 

35.74 

36.77 

0.5429 He - 0.2189 Ne - 0.2382 Ar 

293 

55 and 57 

25.04 

25.90 

0.5000 He - 0.5000 Ar 

5 000 

58 

186.00 

200.30 

0.5000 He - 0.5000 Ar 

10 000 

58 

316.00 

344.20 

0.4000 N - 0.6000 N 2 

10 000 

59 

234.00 

241.80 

0.6000 O - 0.4000 0 2 

10000 

59 

244.00 

240.70 

0.4828 C 2 H 4 - 0.5172 NH3 

523 

56 

17.64 

17.43 

0.6719 H 2 - 0.3281 N 2 

523 

60 

22.02 

23.25 





(a) Comparison with 
reference 45 . 


(b) Comparison with 

references 55 and 56. 


(c) Comparison with 
reference 5^* 


Figure B2.- A comparison of the present proposed relation for mixture viscosity 

with other studies. 


28 



APPENDIX B - Concluded 


When equations (B4), (B5), and (B6) are combined, the mixture viscosity can be 
expressed as 


M = 


(2.2621 X 10-6)t°- 65 
Cj 

Mi^i 




APPENDIX C 


LEWIS NUMBER AUGMENTATION FACTOR 


In the present analysis, the Lewis number augmentation factor 




z i,w) 


is assumed to have a value of 1. This approximation is based on a comparison with the 
equivalent expression for equal diffusion coefficients. The equations for the bifurcation 
approximation of multicomponent diffusion will reduce to the usual equation for binary 
diffusion when the binary diffusion coefficients are assumed equal. 

When the binary diffusion coefficients are assumed equal, then Fj =1, = M > 

and z i = c^, and the approximation for the Lewis number augmentation factor is 



The approximation is exact for a Lewis number of 1 or if the species concentrations do 
not vary through the boundary layer. Also, the approximation is reasonably valid over 
a range of other conditions. 

Let an energy factor E be defined as 


E = 


^^i,w( c i,e _ c i,w) 


hs _ ^w 


or 


E s 


h i,w c i,< 


el ~ ^w 


h s - h w 


Now, 


h i,e = h i,w 


From equations (C3) and (C4), the maximum value of E will be 


(C2) 


(C3) 


(C4) 


®max _ 


hs ~ hw 
h s - h w 


= 1 


(C5) 
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APPENDIX C — Continued 


whereas the minimum value of E will be zero. Therefore, 

OgESl (C6) 

Expressions for the Lewis number augmentation factor are presented in refer- 
ences 8 and 49 for binary diffusion as follows: 


Reference Lewis number augmentation factor 


Present 

i 


1 + 


N Le - 1 „ 

Pl. ) X/ 3 


% 

49 


[1 ♦ (N Le - 1)E]°' 6 


8 1 + [(NLe ) 0 ' 52 - i]e 

The Lewis number augmentation factors are presented in table Cl for a range of Lewis 
numbers and energy factors . As shown by the data, the Lewis number augmentation 

TABLE CL- LEWIS NUMBER AUGMENTATION FACTORS 


E 

N Le 

Lewis number augmentation factor 


Present 

Reference 49 

Reference 8 

0.10 

0.6 

0.953 

0.976 

0.977 


.8 

.978 

.988 

.989 


1.0 

1.000 

1.000 

1.000 


1.2 

1.019 

1.012 

1.010 


1.4 

1.036 

1.024 

1.019 

0.25 

0.6 

0.881 

0.939 

0.942 


.8 

.946 

.970 

.973 


1.0 

1.000 

1.000 

1.000 


1.2 

1.047 

1.030 

1.025 


1.4 

1.089 

1.059 

1.048 

0.50 

0.6 

0.763 

0.875 

0.883 


.8 1 

.892 

.939 

.945 


1.0 

1.000 

1.000 

1.000 


1.2 

1.094 

1.059 

1.049 


1.4 

1.179 

1.116 

1.096 

0.75 

0.6 | 

0.644 

0.807 

0.825 


.8 

.838 

.907 

.918 


1.0 

1.000 

1.000 

1.000 


1.2 

1.141 

1.087 

1.075 


1.4 

1.268 

1.170 

1.143 

1.00 

0.6 

0.526 

0.736 

0.767 


.8 

.784 

.875 

.890 


1.0 

1.000 

1.000 

1.000 


1.2 

1.188 

1.116 

1.099 


1.4 

1.357 

1.224 

1.191 
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APPENDIX C - Concluded 


factor is approximately 1 except for the cases where the Lewis number is significantly 
different from 1 and the energy factor approaches the maximum value. The energy fac- 
tor will approach the minimum for high enthalpies at the outer edge of the boundary layer 
unless the wall enthalpies of the species present at the boundary -layer edge are signifi- 
cantly greater than the wall enthalpies of the species present at the wall. 

From analysis of the present computer results, it was noted that the Lewis number 
at the wall did not vary significantly from 1 for the base gases and gas mixtures which 
are presently considered for planetary entry. Significant variation from 1 was noted for 
the gases in which significant dissociation or reactions occurred at the wall. Thus, it is 
concluded that the approximation of 1 for the Lewis number augmentation factor is rea- 
sonably valid for a wide range of conditions and gases. 
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TABLE I.- MOLECULAR WEIGHTS AND TRANSPORT 


PARAMETERS FOR BASE GASES 


Base gas 

M 

a, A 
(a) 

e/k, K 
(a) 

y 

(b) 

My 

n 2 

28.014 

3.798 

71.4 

0.03654 

1.024 

o 2 

32.000 

3.467 

106.7 

.04129 

1.321 

h 2 

2.016 

2.827 

59.7 

.06775 

.137 

He 

4.003 

2.551 

10.22 

.10845 

.434 

Ne 

20.183 

2.820 

32.8 

.07449 

1.504 

Ar 

39.948 

3.542 

93.3 

.04036 

1.612 

co 2 

44.011 

3.941 

195.2 

.02919 

1.285 

nh 3 

17.031 

2.900 

558.3 

.04605 

.784 

ch 4 

16.043 

3.758 

148.6 

.03345 

.537 


a From reference 50. 
^Calculated from equation (37). 




TABLE H.- HEAT- TRANSFER COEFFICIENTS FROM CORRELATION OF PRESENT COMPUTER RESULTS 


Gas composition 
(composition by mass fraction) 

K, 

Error for correlation, percent 

kg/s -m3/2_ atm 1/2 

Average ] 

Maximum 

n 2 

0.1112 

3.6 

10.0 

°2 

.1201 

3.4 

10.7 


.0395 

3.6 

14.1 

He 

.0797 

6.0 

18.8 

Ne 

.1474 

3.5 

10.5 

Ar 

.1495 

2.7 

7.5 

co 2 

.1210 

3.4 

8.4 

nh 3 

.0990 

9.6 

18.3 

ch 4 

.0807 

8.7 

22.6 

0.23200 2 -0.7680N 2 (air) 

.1113 

3.3 

9.8 

0.6714 CO 2 -0. 3286 Ar 

.1262 

2.6 

6.8 

0.1345 CO 2 -0. 8655 N 2 

.1113 

3.0 

9.4 

0.5000 C0 2 -0.5000 N 2 

.1131 

2.7 

7.4 

0.8500 C0 2 -0.1500 N 2 

.1173 

3.0 

8.5 

0.7500 O 2 -0. 2500 N 2 

.1164 

3.0 

8.8 

0. 5000 O 2 -0. 5000 N 2 

.1132 

3.1 

7.7 

0.1500 H 2 -0.8500 He 

.0657 

4.3 

15.8 

0.3500 H 2 -0.6500 He 

.0547 

3.3 

13.8 

0.6500 H 2 -0.3500 He 

.0459 

3.4 

14.2 

0.5000 H 2 -0.5000 Ar 

.0622 

4.9 

13.1 

0.1500 N 2 -0.8500H 2 

.0473 

5.7 

16.7 

0.5000 N 2 -0.5000 H 2 

.0636 

6.4 

17.8 

0.2000 C0 2 -0.8000H 2 

.0537 

7.6 

23.5 

0.4000 CO 2 -0. 6000 H 2 

.0629 

7.6 

28.9 

0.6000 C0 2 -0.4000H 2 

.0771 

8.6 

27.0 

0.8000 C0 2 -0.2000 H 2 

.1002 

10.2 

17.8 

0.3626 C0 2 -0.3077N 2 -0.3297 Ar 

.1218 

2.7 

9.1 

0.1339 C0 2 -0.8525N 2 -0.0136Ar 

.1118 

2.9 

9.6 

0.3000 Ne-0.3000 Ar-0 .4000 He 

.1143 

6.3 

19.9 

0.3500 H 2 -0.3500 C0 2 -0.3000 N 2 

.0809 

9.0 

21.5 

0.352 H 2 -0.423He-0.176Ne-0.024CH 4 -0.025NH 3 

.0662 

6.2 

16.3 
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TABLE IH.- COMPARISON OF SUMMATION RELATIONS FOR GAS MIXTURES 


Gas composition 

Nr. 


K, kg/s-m^ 

/2- atm 1/2 


(composition by mass fraction) 

lv Pr,w 

(a) 

Computer results 

Present, eq. (43) 

Present, eq. (44) 

Ref. 28, eq. (46) 

n 2 

0.680 

0.1112 




°2 

.685 

.1201 




H 2 

.675 

.0395 




He 

.667 

.0797 




Ne 

.667 

.1474 




Ar 

.667 

.1495 




co 2 

.685 

.1210 




nh 3 

.520 

.0990 




ch 4 

.608 

.0807 




0.2320 O 2 -0.7680N 2 (air) 

.695 

.1113 

0.1117 

0.1131 

0.1131 

0.6714 CO 2 -0. 3286 Ar 

.711 

.1262 

.1252 

.1285 

.1291 

0.1345 CO 2 -0. 8655 N 2 

.692 

.1113 

.1113 

.1124 

.1124 

0.5000 CO 2 -0. 5000 N 2 

.695 

.1131 

.1145 

.1158 

-1159 

0.8500 CO 2 -0. 1500 N 2 

.691 

.1173 

.1187 

.1194 

.1194 

0.7500 O 2 -0. 2500 N 2 

.698 

.1164 

.1162 

.1177 

.1177 

0.5000 0 2 -0.5000 N 2 

.702 

.1132 

.1134 

.1154 

.1155 

0.1500 H 2 -0. 8500 He 

.697 

.0657 

.0644 

.0659 

.0691 

0.3500 H 2 -0. 6500 He 

.710 

.0547 

.0536 

.0553 

.0588 

0.6500 H 2 -0.3500 He 

.705 

.0459 

.0448 

.0460 

.0480 

0.5000 H 2 -0. 5000 Ar 

.496 

.0622 

.0649 

.0540 

.0625 

0.1500 N 2 -0.8500 H 2 

.570 

.0473 

.0469 

.0424 

.0437 

0.5000 N 2 -0.5000H 2 

.477 

.0636 

.0649 

.0526 

.0583 

0.2000 C0 2 -0.8000H 2 

.540 

.0537 

.0498 

.0436 

.0456 

0.4000 C0 2 -0.6000H 2 

.479 

.0629 

.0606 

.0493 

.0541 

0.6000 C0 2 -0.4000H 2 

.439 

.0771 

.0752 

.0580 

.0662 

0.8000 C0 2 -0.2000H 2 

.457 

.1002 

.0937 

.0740 

.0857 

0.3626 C0 2 -0.3077N 2 -0. 3297 Ar 

.710 

.1218 

.1213 

.1246 

.1255 

0.1339 C0 2 -0.8525N 2 -0.0136 Ar 

.693 

.1118 

.1115 

.1127 

.1128 

0.3000 Ne-0.3000Ar-0.4000 He 

.485 

.1143 

.1275 

.1053 

.1104 

0.3500 H 2 -0.3500 C0 2 -0.3000 N 2 

.430 

.0809 

.0795 

.0606 

.0692 

0.352 H 2 -0.423 He-0. 176 Ne-0.024 CH4-O.O25 NH3 

.515 

.0662 

.0669 

.0571 

.0626 


a Average value from computer results for a range of conditions. 
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TABLE IV.- COMPARISON OF SUMMATION RELATION 
WITH REFERENCES 14 AND 15 


Gas composition 

K, kg/s-m 

V2. alm l/2 

(composition by mass fraction) 

Refs. 14 and 15 

Present, eq. (44) 

co 2 

0.1477 


N 2 

.1259 


Ar 

.1596 


0.8655 N 2 - 0.1345 C0 2 

.1278 

0.1283 

0.3889 N 2 -0.6111 C0 2 

.1369 

..1379 

0.6717 CO 2 -0. 3283 Ar 

.1513 

.1513 

0.5614 Ar-0. 3149N 2 -0. 1237 C0 2 

.1459 

.1450 



Wall temperature, K 


300.0 

300.0 

300.0 

300.0 

300.0 

300.0 

1111.0 

1111.0 

1111.0 

1111.0 

1111.0 

1111.0 


h$ - hyY» MJ/Kq 


q^R/Pj, Mw/m 3/2 -atm 1/2 


(a) Key for all parts of figure 1. 


Symbol 

Pressure 

O 

8 

.001 
.010 
• 100 

A 

1.000 

10.000 

K 

100.000 

Q 

.001 

O 

.010 

0 

.100 

h 

1.000 

® 

10.000 


100.000 



h s - h w 

(b) Nitrogen. 

Figure 1.- Heat-transfer results from present computer study. 



(d) Hydrogen. 
Figure 1.- Continued. 
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(g) Argon. 



20 40 60 80 

h s - h w 

(h) Carbon dioxide. 

Figure 1.- Continued. 



(i) Ammonia. 



(j) Methane. 
Figure 1.- Continued. 



(k) 0.2320 02-0.7680 N 2 (air). 
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(Z) 0. 67 14 CO2 -0.3286 Ar. 
Figure 1.- Continued. 
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(m) 0.1345 CO 2-0.8655 N2. 
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(n) 0.50C02-0.50 N 2 . 
Figure 1.- Continued. 




(q) 0.50 02-0.50 N 2 . 
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(r) 0.15 H 2 -0.85 He, 
Figure 1.- Continued. 




(t) 0.65H 2 -0.35He. 
Figure 1.- Continued. 
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(u) 0.50 H 2 -0.50 Ar. 



(v) 0.15 N 2 -0.85 H 2 . 
Figure 1.- Continued. 
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(w) 0.50N2-0.50H 2 . 



(x) 0.20C02-0.80H 2 . 

Figure 1.- Continued. 
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(ee) 0.35H 2 -0.35C02-0.30N 2 . 



(ff) 0.352 H 2 -0.423 He-0. 176 Ne-0.024 CH 4 -0.025 NH 3 
Figure 1.- Concluded. 




(a) Comparison with equation (42). 



(b) Comparison with equation (41). 

Figure 2.- A comparison between the present analytical analysis and the present computer 
results for the heat-transfer coefficient, (units for K are kg/s-m^/^-atm 



Data from literature 



Figure 3.- A comparison of the present analytical analysis with other studies 


P op 






P s , atm p $t atm 


(a) 0.2320 O 2 -0.7680N 2 (air). T w = 1111 K; (b) 0.85CO2-0.15N 2 . T w = 300 K; 

h s = 81.354 MJ/kg. h s = 104.596 MJ/kg. 

Figure 5.- The trend of the heat-transfer coefficient with pressure, (units for K are kg/s-m 3 / 2 -atm 1 / 2 . 
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(a) 0.2320 0 2 -0.7680N2 (air). T w = 300 K; 
p g = 1.0 atm. 




h s , Mi/kg 

(d) 0.3626 CO 2 -0.3077N2-0.3297Ar. 
T w = 300 K; p s = 1.0 atm. 
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(e) 0.65 H 2 -0.35 He. T w = 300 K; 
p s = 0.001 atm. 



h s , MJ/kg 

(c) Carbon dioxide. T w = 1111 K; p s = 0.1 atm. 


.14 
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h s , MJ/kg 

(f) 0.6714 CO2-0.3286Ar. T w = 1111 K; 
p s = 10.0 atm. 


Figure 6.- The trend of the heat-transfer coefficient with enthalpy. 


(units for K are kg/s- mV2_atml/2.) 
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(a) CO2-N2 mixtures. (b) He-H 2 mixtures. 
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(c) CO2-H2 mixtures. (d) N2-H2 mixtures. 

Figure 7 .- A comparison of the summation relations with the present computer results 
for some binary gas mixtures, (units for K are kg/s-m 3 / 2 -atml/ 2 .) 
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M 0 

Figure 8.- A comparison of the present computer results with the molecular weight 
correlation of reference 16. (units for K are kg/s -mV2- atm 1/2 
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